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Abstract 

A numerical investigation is reported of turbulent incompressible jets confined in 
two ducts, one cylindrical and the other conical with a 5° divergence. In each case, 
three Craya-Curtet numbers are considered which correspond, respectively, to flow 
situations with no, moderate and strong recirculation. Turbulence closure is achieved 
by using the K-e model and a recently proposed realizable Reynolds stress algebraic 
equation model that relates the Reynolds stresses explicitly to the quadratic terms 
of the mean velocity gradients and ensures the positiveness of each component of the 
turbulent kinetic energy. Calculations are carried out with a finite- volume procedure 
using boundary-fitted curvilinear co-ordinates. A second-order accurate, bounded 
convection scheme and sufficiently fine grids are used to prevent the solutions from 
being contaminated by numerical diffusion. The calculated results are compared 
extensively with the available experimental data. It is shown that the numerical 
methods presented are capable of capturing the essential flow features observed in 
the experiments and that the realizable Reynolds stress algebraic equation model 
performs much better than the K-c model for this class of flows of great practical 
importance. 


1 Introduction 

This paper is concerned with numerical computations of confined turbulent jets. 
The general features of the flow under consideration are sketched in Fig.l. At 
the entrance, two uniform flows, a jet of high velocity and an ambient stream of 
low velocity, are discharged into the duct. Due to turbulent entrainment, the jet 
increases its mass flux while spreading and this causes an equal decrease in the 
mass flux of the ambient stream. An adverse pressure gradient is thus set up by 
the decrease in the velocity of the ambient stream which can be considered as a 


1 



potential flow. When the ratio of jet to ambient velocities at the entrance is above 
a critical value, a recirculation zone occurs at the duct wall downstream of the inlet 
plane. This is because the jet has consumed the whole ambient flow before reaching 
the wall and further entrainment must create reverse flows in order to maintain the 
total mass flux conservation. For a given geometry, recirculation as well as adverse 
pressure gradients can be intensified by increasing the jet to ambient velocity ratio 
at the entrance. These features are common to many flows involving two streams of 
different velocities in industrial apparatus, in particular, in combustion chambers in 
which recirculating currents are used to enhance flame mixing and stability, and in 
ejectors whose role as a jet pump is achieved through the pressure rise generated by 
the jet entrainment. Therefore, the understanding of confined jet flows is of great 
practical importance to the design of such devices. 

A large amount of theoretical and experimental research has been devoted to 
confined jets (Rajaratnam, 1976). Early theoretical studies were restricted to the 
development of integral methods based on similarity assumptions ( Cray a and Curtet, 
1955, Curtet, 1960, Hill, 1967 and 1973, Mikhail 1960 and Razinsky and Brighton, 
1972). The essential features of these methods are presented in detail in Rajaratnam 
(1976). Among the most well-known of them is the approximate confined jet theory 
of Craya and Curtet. They found that confined jets in constant area ducts can form 
a family characterized by a single non-dimensional parameter C t termed the Craya- 
Curtet number, provided that the flows remain self-similar. This condition cannot 
be satisfied in general, especially in the presence of recirculation, and the Craya- 
Curtet number is not constant in variable area ducts. Nevertheness, the values of C t 
at the entrance remain useful to characterize the inlet flow conditions. Experimental 
studies can be found in Barchilon and Curtet (1964), Barchilon (1969), Binder and 
Kian (1983), Kian (1981) and Razinsky and Brighton (1971). Experiments show 
that recirculation occurs in cylindrical ducts when C t < 0.96 and in a conical duct 
with a 5° divergence when C t < 1.1. 

The general applicability of the classical integral method is invalidated by the 
complexities of the flow. At the present time, the only feasible theoretical solution 
is to use the differential method that solves the complete governing equations with 
the aid of turbulence modeling. Several authors (Gosman et al., 1979, Habib and 
Whitelaw, 1979 and 1981, Jones and Marquis, 1985, Khalil et al., 1975 and Zhu, 
1986) have reported the relevant calculations using such differential methods. In 
these calculations, turbulence effects were represented either by the K-e model or 
by second-moment closures. The success of the differential methods depends, to 
a large extent, on the performance of the turbulence model used. These numeri- 
cal studies showed that the calculations with the K-e model could reproduce some 
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basic flow characteristics observed in the experiments in general, but large discrep- 
ances were present in detail. These discrepances were mostly attributed to inherent 
imperfections of the model such as the isotropic eddy-viscosity representation and 
the insensitivity to streamline curvature effects. Second-moment closures have no 
such drawbacks and represent a higher level of turbulence modeling beyond the 
eddy-viscosity level. However, it was found in a systematical numerical study (Zhu, 
1986) that second-moment closures, whether of the algebraic or of the transport 
form, did not manifest themselves to be significantly superior to the K-e model in 
the computation of the confined jets. It should be pointed out that the previous 
calculations were almost exclusively performed on coarse grids and with the hybrid 
(central/upwind differencing) scheme (Spalding, 1972) that is essentially of first- 
order accuracy in convection-dominated situations. Recent work of Zhu and Rodi 
(1992) has shown that the calculations of confined jets were very sensitive to the 
approximation of the convection terms of transport equations and the hybrid scheme 
could not yield grid-independent results, especially for the turbulent quantities, on 
a grid as fine as 102x82 points. Therefore, it is certain that the early calculations 
were contaminated to some extent by numerical diffusion and the performance of the 
turbulence models for the confined jet flows needs to be re-assessed on a numerically 

reliable basis. 

Recently, Shih and Lumley (1993) have shown that the Reynolds stress, being 
a second rank tensor, can be expressed as a fourth-order polynomial of the mean 
velocity gradients. This is the most general stress-strain relationship within the 
framework of algebraic turbulence modeling, with the linear stress-strain relation in 
the Boussinesq’s eddy-viscosity concept being its first-order approximation. Based 
on this, Shih et al. (1993) have recently proposed a quadratic stress-strain relation in 
conjunction with the two modeled equations of K and £ for practical calculations. 
Realizability constraints (Lumley, 1978 and Schumann, 1977) have been used to 
derive appropriate functional expressions for the model coefficients, in particular, 
the coefficient C ^ has been naturally related to the time scale ratio of the turbulence 
to the mean strain rate. As a result, the model ensures the positivity of individual 
turbulent normal stresses, an important feature that is not present in most existing 
turbulence models. The empirical constants in the model have been fine-tuned in the 
calculation of two backward-facing step flows. As the model appears promising as a 
competitive alternative in the turbulence modeling arsenal, it is of interest to test its 
performance for confined jets which, due to their complicated flow structures such 
as strong recirculation coupled with severe adverse pressure gradients, constitute a 
challenge to the turbulence modeling. 

To this end, calculations are performed for the flows studied experimentally by 
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Barchilon and Curtet (1964) in a cylindrical duct and by Binder and Kian (1983) 
in a conical duct with 5° divergence. These two experiments have been chosen in 
this study because of their diverse flow conditions which vary gradually from no 
to strong recirculation. The standard K-e model is also used, since it is the most 
popular turbulence model used today in calculations of complicated flows and can 
also be used to highlight the performance of the new model. The numerical accu- 
racy of solutions is ensured by using a second-order accurate, bounded convection 
scheme (Zhu, 1991a) and sufficiently fine grids. The performances of the models are 
examined through extensive comparisons with experimental data. 


2 Calculation Approach 

2.1 Turbulence Models. Incompressible turbulent flows are governed by the 


statistically averaged continuity and Navier-Stokes equations. The Reynolds or tur- 
bulent stresses ( = — puiUj ) appearing in these equations are calculated by using 
the following two models: 

1) K-e model (Launder and Spalding, 1974) 


2 

T ij = T Uj,i) ~ 

(1) 

pK 2 

pt — Cfi y CV — 0.09 

£ 

(2) 

where the turbulent kinetic energy K and its dissipation e are 
tively, by the modeled transport equations given in Table 1. 

calculated, respec- 

2) Realizable Reynolds stress algebraic equation (RRSAE) model (Shih et al., 1993) 

T ij = PtiPij + U j,i ) + Tij — -pKSij 

(3) 

u C PX 1 c 2/3 

(4) 


where 


+CV 2 (tfi.*tfi.fc - + c r3 (u Ki u kJ - (5) 
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n = UkjtUu . , 


n = Uk,iUk,i 


( 6 ) 


, = — , 5 = (25«S ji )'/> , S* = hVij + £/,,) (7) 

e z 

f = ^, n = (MW 1 ". n; ; = W., - f/ ; .i)/2 + 4 £mji u, m (8) 

a> m is the rotation rate of the reference frame, and the model constants are 

CVi = —4, (7.2 = 13, a r3 = -2, A 2 = 1000. (9) 

In the work of Shih et al. (1993) and Zhu and Shih (1993), the following two sets 
of values for A\ , 7 have been tested 

At = 5.5, 7 = 0 (10) 

Ai = 1.25, 7 = 0.9 (11) 

and both were found to give almost identical predictions for the two backward-facing 
step flows. With Eq.( 11), the rotational effect of the mean flow explicitly enters 
into Cp. However, the values in Eq.( 10) are taken in the present work. The K and 
e in the RRSAE model are calculated with the same equations as in the K-e model. 
It can be seen that the RRSAE model reduces to the K-e model if the quadratic 
terms T tJ * are set to zero. 


2.2 Curvilinear Form of Governing Equations. Boundary-fitted curvilinear 
co-ordinates are used in order to handle irregular flow boundaries. In such co- 
ordinates, the two-dimensional, steady-state governing equations using Cartesian 
velocity components Ui can be written in the following general form. 

(Ci<t>-D i4> ), Xi =r a JSt, i= 1,2 ( 12 ) 

where the convective coefficients ( 7 ,, the diffusion terms Di^> y and the source terms 
Sffr are given in Table 1 for different dependent variables <j>\ J is the Jacobian of the 
transformation between the curvilinear coordinates and the reference Cartesian 
coordinates y { . Eq.( 12) includes both the plane (a = 0) and axisymmetric (a = 1) 
forms. In the latter form, y\ (= z) and 3/2 ( — f*) nre in the axial and radial directions, 
respectively. 

In the curvilinear transformation, Eqs. ( 1 ) - ( 8 ) all remain the same, while the 
velocity gradients are calculated as follows: 
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( 13 ) 


U \.i=(P{Vi, Il +P 2 l U l „)/J 
Ul'2 = (/Jjf/i,*, + PlUi'Z,)/ J 
U 2 ,i = {P\U 2l , x +fi(J itX2 )/J 
u 2,2 = (PlU 2 , Xl +(3 2 2 U 2 „)/J 
U 3i3 = aU 2 /r 

where Uj iXi is the partial derivative of Uj with respect to x 

Table 1. Individual terms in the transport equations 
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Ci = pr a {(3\U\ T P\ U 2 ) , C 2 = pr a ((5\{Ji + P 2 Ut) 
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(3j = cofactor of yj tXi in J 

B) = p\P{ + p&i , u>j = P}u im + p*Ui„ , 

2 

P — -| — pK (p is the mean static pressure) 

3 

r n = p+ti } (n = 1,2, 3,4) 

G = TiiUl'l + T 22 U 2i2 4" 0^33^3,3 4 T f^2,l) 

(7 lt = 1.44, (72c = 1.92, <7 j fr 2 -3- <T 2 1, <r 4 ~ 1.3 


7 = 


6 







2.3 Boundary Conditions. Computational boundaries involved are inlet, out- 
let, axis of symmetry and solid wall. Among them, the inlet conditions require 
special attention because they have a considerable influence on the calculations 
(Zhu et al., 1987 and Zhu and Rodi, 1992). At the inlet, the jet and the ambient 
flow have uniform velocities Uj and U a (Fig.l), but the flow conditions in the initial 
shear layer between the jet and the ambient flow are unknown and determination of 
these conditions is not a trivial matter. In the present work, the parabolic entrance 
region (PER) scheme of Zhu et al. (1987) is used. The PER scheme was developed 
on the assumption that although the flow as a whole is elliptic, there is a short 
region near the entrance where the flow is parabolic. A parabolic calculation is first 
carried out over a short distance between x = 0 and x = z c by using the following 
mixing length model 


P UV=-fi t U, r , fi t — C ,2 (l '2 — ri) 2 \U,r | (14) 

where U(= U\) is the axial mean velocity, r x and r 2 are the coordinates of the inner 
and outer edge of the initial shear layer (Fig.l) and C is an empirical coefficient 
derived from the experiment of Rajaratnam and Pani (1972) 

C 2 = 0.0042 + 0.004 UJUj , 0 < U a /Uj < 0.2 (15) 

The results of the parabolic calculation are then used as the inlet conditions at 
x = x e for the elliptic calculation. At this location, the inlet values of K and e are 
calculated by 


K = -tZv/0.3 , e = 0.09 pK 2 /fi t (16) 

It was found (Zhu, 1986) that the PER scheme gives satisfactory predictions in 
the parabolic entrance region and that the elliptic calculations are insensitive to x € 
provided that 1 < x e jd 0 < 3. 

The outlet boundary is placed at x = 10Z? a which is sufficiently far away from 
the main region of interest. At this boundary stream wise gradients of all variables 
are set to zero. Along the axis of symmetry the normal velocity component and the 
normal gradients of the other variables are set to zero. The standard wall-function 
approach (Launder and Spalding, 1974) is used to bridge the viscous sublayer at the 
duct wall. 

2.4 Numerical Procedure. The system of transport equations ( 12) is solved 
by using a conservative finite-volume method designed for calculating incompressible 
elliptic flows with complex boundaries. The method uses a non-staggered variable 
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arrangement, namely, all the dependent variables are stored at the geometric center 
of each control volume. The momentum interpolation of Rhie and Chow (1983) 
is used to avoid checkerboard oscillations usually associated with non-staggered 
grids. The velocity-pressure coupling is achieved via the SIMPLEC algorithm (Van 
Doormal and Raithby, 1984). Numerical implementation of the RRSAE model is 
straightforward; the linear part is treated in the same manner as in the K-e model 
and the quadratic part is included in the source terms. To ensure numerical accuracy 
and stability, the convection terms of all the transport equations are differenced by 
the hybrid linear/parabolic approximation (HLPA) of second-order accuracy (Zhu, 
1991a), and all the other terms by the conventional central differencing scheme. The 
resulting set of algebraic difference equations is solved with the strongly implicit 
procedure of Stone (1968). The iterative solution process is regarded as converged 
when the maximum normalized residue of all the dependent variables is below 0.5%. 
The details of the present numerical procedure are given in Rodi et al. (1989) and 
Zhu (1991b). 

3 Results 

The experiments of Barchilon and Curtet (1964) and Binder and Kian (1983) 
were taken as the test cases, hereinafter termed BC- and BK-case, respectively. The 
geometrical dimensions and the inlet flow conditions of both cases are given in Tables 
2 and 3. 


Table 2. Geometrical dimensions of the ducts 


case 

D 0 (cm) 

d 0 (cm) 

6 (degree) 

L (cm) 

BC 

16 

1.2 

0 


BK 

16 

1.6 

2.5 

64 


In the BK-case, the velocity field and the pressure field were measured separately 
using the different inlet velocities as shown in Table 3. The Craya-Curtet number 
C t is defined by 



where m is the total momentum flux non-dimensionalized with the total flow rate 
Q and the duct area S. 
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Table 3. Inlet flow conditions 


case 

c t 

Uj (cm/s) 

U a (cm/s) 


0.976 

1293.6 

84.48 

BC 

0.506 

1253.8 

39.81 


0.152 

1296.2 

7.42 

BK 

1.230 

40 

4.885 

(velocity 

field) 

0.775 

40 

2.662 

0.590 

40 

2.330 

BK 

1.090 

650 

61.38 

(pressure 

field) 

0.842 

650 

46.58 

0.570 

650 

29.90 


Calculations were made on tne oray uvn — ” T , ttv 

the solutions was examined by using the two convection schemes, HLPA and HY- 
BRID (central/upwind differencing), and three grids consisting of 50x40 (grid 1), 
86x50 (grid 2) and 120x80 (grid 3) points, respectively. The HYBRID scheme tha 
is highly diffusive in the presence of both convective dominance and flow-to-gnd 
skewness has been used here only to highlight the importance of using ^gher-ordcr 
accurate schemes. Test results obtained with the RRSAE model^at C, - ■ 

in the BC-case are shown in Fig.2 for the turbulent shear stress nu-proRles at the 
downstream location */D.= 1.875. It can be seen that the results with HLPA on 
the coarse grid 1 are already very close to those on the fine grid 13,' while significant 
differences exist between the corresponding results with HYBRID. The HLPA re- 
sults on the intermediate grid 2 can be considered as grid-independent because the 
refinement from the grid 2 to the grid 3 produced differences too smaU to be seen 
on the graph The HYBRID solutions, however, responded to the grid refinemen 
in such a slow manner that they stiff have not reached the grid-independent stage 
on the finest grid. The HYBRID scheme is, therefore, inappropriate for compute- 
rs of confined jets in presence of recirculation. The iterations and CPD-time in 
minutes required for the calculations with HLPA were 196 and 0.2 on grid 1, 640 
and 1.4 on grid 2 and 1874 and 9.3 on grid 3. The calculations with HYBRID took 
about 0.6-0.8 of these numbers. The grid 2 and HLPA were used for all subsequent 

calculations. 

Figs 3 and 4 show the decay of the centerline velocity £/„, normalized by the 
mean velocity of the section U m , at different C, numbers. The figures clearly reveal 
the existence of the potential core characterized by the constant U a . Beyond the 
potential core, U a decays quickly, especially at small C t numbers. The experiment 
data show that the potential core is shorter in the BK-case than in the BC-case, 
because the former has a larger mixing layer between the jet and the ambient flow 
at the entrance. This feature is also reflected in the computed results. The RRSAE 
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model predicts a longer potential core length than the K-e model does, but it is 
difficult to judge which is better because this length cannot be precisely determined 
from the first and second experimental points. Overall, both models predict the 
centerline velocity decay well, though small underpredictions occur in the K-e model 
results at the intermediate and small values of C t in both cases. Figs. 5 and 6 show 
the axial mean velocity profiles at four downstream locations, at the smallest value 
of Ct in each case. In the BC-case, the RRSAE model results are in excellent 
agreement with the experimental data whereas noticeable discrepancies exist in the 
K-e model results at the last three downstream locations. In the BK-case, the results 
of RRSAE model are also better than those of the K-e model, but the agreement is 
less satisfactory than in the BC-case. At the location x/D 0 = 1.25, the experimental 
profile shows no constant ambient velocity portion while the calculated profiles by 
both models still have a vertical plateau around rj R — 0.6. The comparison at the 
location xjD 0 = 2.5 indicates that the width of the predicted reverse-flow region is 
somewhat too thin. Comparing the corresponding profiles in the BC-case (Fig. 5) 
which has a much stronger recirculation, the experimental velocity minimum at 
the last location in Fig.6 seems to be too far away from the duct wall. It is to 
be noted that the flow is highly perturbed in the recirculation region in which the 
measurement uncertainty is likely to be greatest. 

In order to study the jet spreading, the excess flow rate Qj has been introduced 
which is defined as follows 


Qj = 2tt [ T, (U — Ui)rdr (18) 

Jo 

where U\ is the ambient velocity (Fig.l). In the recirculation region, the ambient 
velocity has no physical meaning and is defined as the minimum velocity (having 
negative value) for analytical convenience. Figs. 7 and 8 show the variation of the 
excess flow rate Qj with x and Ct. As a consequence of the turbulent entrainment, 
the excess flow rate increases first, passes through a maximum and then decreases 
in case of recirculation which corresponds to Qj/Q > 1 . This variation becomes 
more pronounced as Ct decreases. The calculations agree well with the experiments 
at large C t , but the agreement deteriorates as recirculation intensifies. At small C t 
(<0.6), the calculations overpredict the maximum excess flow rate in the BC-case, 
while the opposite happens in the BK-case. It should be pointed out that the excess 
flow rate, due to its definition, is a quantity that is highly sensitive to errors in the 
velocity profiles so that a small change in U\, especially in the recirculation zone, 
will result in a large difference in Qj. Regarding the comparison between the two 
models, the RRSAE model clearly performs better than does the K-e model. 

The separation and reattachment points are given in Table 4. They arc the 
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locations where Q j j Q—\ in Figs. 7 and 8. In the BK-case, only a range of 3.4 ~ 3.8 D a 
was given experimentally for the reattachment point due to the high unsteadiness of 
the recirculating bubbles. The experiments in both the BC- and BK-cases revealed 
that as C t decreased, the separation points moved upstream while the reattchment 
points remained basically unchanged. The RRSAE model captures this feature well 
and predicts the locations of the recirculation bubbles better than the K-e model. 


Table 4. Separation and reattachment points ( x,/D 0 , x r /D 0 ) 


case 

c t 

experiment 

K-e model 

RRSAE model 

BC 

0.506 

1.70, 3.07 

1.37, 2.7 

1.6 , 3.03 


0.152 

0.45, 3.07 

0.40, 2.95 

0.5 , 3.15 

BK 

0.775 

2.50, 3.4~3.8 

1.82, 3.17 

2.15, 3.79 


0.590 

1.50, 3.4~3.8 

1.43, 3.22 

1.45, 3.81 


Figs.9 and 10 show the variation of pressure coefficients along the duct walls. 
Here, C p is defined by 


, A p - pUlj2 

pUj/% 


(19) 


and A p is the pressure difference between the location x and the entrance, dhe 
pressure gradient is governed by the jet entrainment, the contraction and expansion 
of the flow caused by recirculating eddies as well as the geometry of the duct. 
As shown in Kian (1981), the entrainment and the divergence of the duct can only 
produce a maximum pressure difference equal to pU\l 2, while the pressure difference 
created by the divergence of streamlines in the downstream part of the recirculating 
bubble is of the order of 


Ap 

Wi 


(1 + K 


>(£)' 


( 20 ) 


where the coefficient K c is a positive quantity defined by the flow velocity, and R c 
is the minimum radius of the contracted stream. Therefore, the pressure rise in the 
recirculating zone depends on the width of the recirculating bubble and can be much 
larger than pU^/2. This explains the variation of C p with C t seen in Figs.9 and 10. 


Regarding the comparison between predictions and experiments, it can be seen that 
although both models predict practically the same total pressure rise which is in 
very good agreement with the measurements (Fig. 9), the RRSAE model captures 
the location where the pressure starts to shoot up much better than the K-e model 
at all the G t numbers in both the cases. This location corresponds to the abscissa 
of the maximum excess flow rate shown in Figs. 7 and 8. 
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Detailed experimental data for the turbulent stresses tZu, WO and uv are available 
only in the BK-case at Ct — 0.59. The computed and measured radial profiles of 
these quantities at four downstream locations are compared in Figs. 11-13. With 
regard to the turbulent normal stresses shown in Figs. 11 and 12, the experimental 
data are basically followed by the results of both the RRSAE model and the K-e 
model, with the former predicting more anisotropy than the latter. The experimen- 
tal data at x/D 0 — 2.5 are seen to exhibit a different trend for both ITu and vv, 
which may possibly be due to measurement errors. The flow visualization in the ex- 
periment indicated that the global flow pattern was highly unsteady in the presence 
of recirculation. With due regard to flow complexities and measurement difficulties, 
the agreement between the predictions and measurements seen in Figs. 11 and 12 
should be considered as reasonably good, but it is difficult to judge which model 
performs better for the turbulent normal stresses, overall. For the turbulent shear 
stress uv shown in Fig. 13, the results obtained with the RRSAE model are clearly 
better than those with the K-e model for all the locations considered. The large dis- 
crepancy seen at x / D a = 2.5 is partially due to the underprediction of the width of 
the backflow region and partially due to the experimental uncertainty, as evidenced 
by the fact that in the experimental data, the change in sign of the shear stress 
profile occurs much further away from the duct wall than the velocity minimum. 


4 Concluding Remarks 

A numerical study has been conducted on axisymmetric confined jets in two ducts 
with and without divergence. The focus of the study was on the performance of the 
two turbulence models under various flow conditions ranging from no to strong 
recirculation. The RRSAE model is new, and its capability to predict this type of 
flows needs to be tested. The K-e model was applied to the same problem before, but 
we believe that a re-assessment is necessary because the previous studies were made 
on an inaccurate numerical basis (first-order differencing together with very coarse 
grids). This study has shown that the confined jet calculations are very sensitive 
to the convective approximation of the transport equations. In order to achieve a 
numerically credible solution without involving excessive grid points, higher-order 
differencing schemes have to be adopted. In this regard, the second-order accurate 
HLPA scheme behaves quite well and its solution is virtually the same as those of 
third-order accurate SMART and SHARP schemes (Zhu, 1992), all three being of 
oscillation-free nature. Therefore, it can be concluded with a sufficient degree of 
confidence that the present calculations reflect the real predictive capability of the 
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models for the confined jets. 

The calculations have been extensively compared with the experiments. The 
comparison clearly shows the superiority of the RRSAE model over the K-c model 
under all the circumstances considered. Specifically, the former performs much bet- 
ter than the latter in the predictions of the pressure distribution, the separation and 
reattachment points of recirculation region, the excess flow rate, the jet expansion 
and the turbulent shear stress. Only for the maximum reverse flow rate, an impor- 
tant parameter to characterize the performance of combustion chambers, does the 
RRSAE model result in little improvement. Since this quantity is directly related to 
the flow near the wall in the recirculation region, the use of the wall function may 
constitute a source of error and also measurement errors should not be excluded in 
this highly perturbed region. Overall, the RRSAE model is shown to be able to 
predict confined jets with an accuracy sufficient for engineering applications. 
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Figure 1. Flow configuration and notations 
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Figure 2. Numerical accuracy test at C t — 0.506 in the BC-case 
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Figure 3. Centerline velocity decay in the BC-case 
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Figure 5. Axial mean velocity profiles at Gt — 0.152 in the BC-case 

(Notation as in figure 3) 
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Figure 6. Axial mean velocity profiles at C t = 0.59 in the BK-case 
(Notation as in figure 3) 







Figure 9. Pressure coefficient along the duct wall in the BC-case 
^Notation as in fieure 31 
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Figure 11. Turbulent normal stress tru- profiles at C t = 0. 

(Notation as in figure 3) 

59 in the BK-case 
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